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In a recent paper L.G. Moretto et al Q claim that the negative heat capacities presented in our 
previously published paper are "artifacts" coming from the use of periodic boundary conditions 
in the Lattice-Gas calculations. We stress in this comment that this claim is wrong: in ref. we 
did not use periodic boundary conditions and anyhow the boundary conditions are irrelevant for the 
statistical ensemble used in Q. The second claim of is that, because of the Coulomb repulsion, 
systems "with A > 60 should present no anomalous negative heat capacities". We show that this 
conclusion is contradicted by exact Lattice-gas simulations including Coulomb forces which present 
negative heat capacities even for A > 200. 



Let us start with the discussion about the boundary 
conditions. It is clearly stated in ref. 2} thaf'L 3 [the lat- 
tice size] is large enough (typically greater than 20 3 lat- 
tice sites) so that the boundary conditions do not affect 
the calculations with a constraining A." This means that 
the results of ref. do not depend upon the conditions 
used at the boundary, and in fact the calculations were 
made without periodic boundary conditions. At that time 
we checked the independence of the boundaries compar- 
ing N = 8000 and N = 27000 lattices. To see a sizeable 
effect of the boundary, we have decreased the size of the 
lattice. A fast calculation of 50000 events for 216 parti- 
cles in a N — 5832 lattice at an energy E — 0.4e confined 
by a Lagrange multiplier log A = —8 , gives a temperature 
T = (0.687 ± 0.004)e, a heat capacity C = -16.3 ± 0.2, 
and a kinetic energy fluctuation cr\/T 2 = 1.64 ± 0.03. If 
periodic boundary conditions are imposed, the temper- 
ature becomes T = (0.682 ± 0.004)e, the heat capacity 
C = -16.9±0.4 and the fluctuation a 2 k /T 2 = 1.63±0.07. 
The temperature decrease in the phase transition region 
is AT = (7.6 ± 0.4)10~ 3 e (AT = (6.7 ± 0.6)10~ 3 e) with 
(without) periodic boundary conditions. This shows that 
the boundary conditions do not affect the thermodynam- 
ics. Indeed, in ref. Q we have analyzed an ensemble 
of particles in the vacuum for which only the average 
spatial extension of the system is defined (isobar ensem- 
ble). Then, provided that the lattice in which the sys- 
tem has still to be discretized for technical reasons is big 
enough, boundary conditions are irrelevant since the par- 
ticles never explore the outer region. 

In a previous paper || we have considered a canonical 
Lattice Gas model in a box of constant volume (isochorc 
ensemble), and there we have used periodic boundary 
conditions. In these calculations negative compressibil- 
ities are reported but the heat capacity is always posi- 
tive. In fact, in the canonical ensemble the heat capacity 
is proportional to the energy variance and can never be 
negative p|. 

Concerning the isochore microcanonical ensemble, the 
published results |5j show that the heat capacity at con- 



stant volume Cy is always positive. Therefore, contrary 
to what is claimed in ref. published results dealing 
with Lattice Gas calculations with typical nuclear size 
systems and periodic boundary conditions 0, @ do not 
show any negative heat capacity. 

We may incidentally note that recently a debate has 
started in the statistical physics community after the ob- 
servation of a negative heat capacity branch for very large 
systems in the the constant-magnetization microcanoni- 
cal Ising model Q (equivalent to the microcanonical iso- 
chore Lattice Gas model). In these cases however, the 
densities involved are so low that once again boundary 
conditions are irrelevant. Moreover, the sizes discussed 
are orders of magnitude larger than the nuclear ones. 

Finally we would like to stress that we have recently 
demonstrated that negative heat capacities in finite sys- 
tems are the origin of first order phase transitions with a 
finite latent heat: to present an energy discontinuity at 
the thermodynamic limit, finite systems must present a 
negative microcanonical heat capacity if the number of 
particles is large enough yj. 

Let us now comment about the effect of Coulomb. The 
conclusion of ref.[l| is based on a simple model of a unique 
cluster in equilibrium with a gas. However, this config- 
uration is not the most probable one both in statisti- 
cal models [jj 0] and in nuclear physics experiments at 
high excitation energy ^(j- I n particular, the negative 
heat capacity region is characterized by a large (m > 3) 
multiplicity of fragments |llj |. Therefore any conclusion 
drawn considering only the specific channel of multiplic- 
ity 1 is not valid. If we would follow the simple model of 
ref. [lj and consider partitions with m = 3 — 5 fragments 
of similar mass [13, this would lead to a figure analogous 
to figure 5 of ref. [l| but with an abscissa approximately 
multiplied by m (or more, since all the light particles have 
to be included), because the inter- fragment coulomb in- 
teraction can be neglected if the volume is large. One 
would then find A « 200 as the limit of negative heat 
capacities. 

Of course fragments do not have the same size and 
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FIG. 1: Caloric curve for the Lattice Gas hamiltonian aug- 
mented with the Coulomb interaction in the microcanoni- 
cal isobar ensemble (log A = -7, A = 207, Z = 82) from 
ref.|l4|. The coupling between neighboring sites is fixed to 
e = 5.5MeV. 



to get a quantitative result one needs to compute the 
relative weight of all the possible channels: this is indeed 
what is done in nuclear statistical models which 
show that the intermediate state between the compound 
nucleus and the vaporized system is multifragmentation: 
the presence of many drops in equilibrium. 

In the case of the Lattice Gas model the different parti- 
tions and their relative weight can be calculated without 



any approximation. In the negative heat capacity energy 
regime there are in average three fragments of size greater 
than 4 2]. The sudden opening of this multifragment 
channel|9j causes a convex intruder in the entropy that 
is responsible of the negative heat capacity. The effect of 
the Coulomb interaction in the Lattice Gas model can be 
appreciated from oFigure 1 which shows a caloric curve 
at constant A = Bp for a system of A = 207 particles and 
charge Z = 82 14]. A clear backbending is visible. In 
the charged Lattice-gas model, the energy interval cor- 
responding to the backbending decreases with increasing 
charge [l4(, showing that the idea that Coulomb tends 
to suppress the negative heat capacity is qualitatively 
correct. However, when all the available channels are 
correctly weighted, the quantitative effect is very differ- 
ent then the one reported in ref. |lj . Only a slight effect 
of the Coulomb interaction has also been reported in the 
isochore framework 

Obviously the result of Figure 1 is model dependent. 
Different macroscopic models show a higher [TsT ] or lower 
[r2l | sensitivity to the Coulomb depending on the de- 
tailed implementation of the Coulomb and nuclear inter- 
actions. On the other hand molecular dynamics (which 
can be solved exactly as the Lattice Gas model) leads to 
caloric curves which are almost independent of Coulomb 
[l6[ . Whether this can be interpreted as an evidence of 
metastable long time tails, the relaxation time of long 
range interactions being excessively long, is a subject of 
debate 01- I n an y case these examples and the results 
of Figure 1 show that a definitive understanding of the 
effect of a non saturating long range interaction on a 
first order phase transition cannot certainly be achieved 
through the oversimplified model of ref. 
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